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This  paper  reviews  known  limit  theorems  for  the  population 
siaes  of  a Bienaymtf«Galton-Watson  process  allowing  Immigration. 

For  the  non* critical  cases  It  is  known  that  the  limit  distribution  is 
rjon-defective  iff  a logarithmic  moment  of  the  Immigration  distribution 
is  finite.  The  new  results  of  this  paper  are  concerned  with  the  sltuatiott 
where  this  moment  is  Infinite  and  give  limit  theorems  for  a certain 
slowly  varying  function  of  the  population  slae.  A parallel  diacusslon  la 
given  for  the  critical  case  and  also  for  the  continuous  time  process. 

The  methods  of  the  paper  are  used  to  give  some  results  on  the 
rate  of  decay  of  the  transition  probabilities  and  on  the  growth  rate  of 
the  stationary  measure.  These  in  turn  are  used  to  obtain  some  limit 
theorems  for  a reversed  time  process.  Known  results  and  applications 
are  reviewed. 


Keywords:  Blenaymd-Galton- Watson,  ixxunlgration,  limit 

theorem,  non-linear  norming,  hlarkov  branching 
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1 . Introduction 

We  ehnil  consider  the  Bienaym4>Galton« Watson  process  allowing 
immigration  (BPI)  which  is  a Markov  chain,  nsO,l, . . . with 

minimal  state  space  j c ^ {0, 1 , . . . ] and  one  step  transition  probabilities 
given  by 


p^^  * coeff.  of  s^  in  h(sHf(*))^ 


(i.JcK) 


where  h(s)aS^i®^  •“d  £(s)«2p.»^  (S(»).  ® S (•)!)  are  probability  gen- 
J J ^ j > 0 •’ 

e rating  functions  (pgf's)  and  h^<l,  p ^ / 1 • The  process  x i**  usually 
thought  as  describing  the  evolution  of  the  size  of  a population  in  which  each 
Individual  produces  j progeny  with  probability  p^  at  (or  by)  the  end  of  its 
life  and  time  is  measured  in  generations.  With  probability  h^,  j immigrants 
enter  the  nth  generation  azid  contribute  to  the  (n-(’l)th  generation  in  the  same 
way  as  natives.  All  individuals  reproduce  independently  of  each  other  and  of 
the  inunigration  process,  and  the  numbers  immigrating  into  successive  gen- 
erations are  independent.  Thus  X represents  the  total  number  of  individuals 

n 

in  the  nth  generation.  This  intuitive  picture  Implies  the  recursive  relationship 


r . ^ . I 

n+1  ^i  ^n+1 


(n«K) 


(1) 


where  * i t lN^,n«  N}  (W^  = [1,2,...])  are  independent  and  have  pgf 

f(*)t  (I^tnelN^}  are  independent  with  pgf  h(*)  and  both  seijuences  are 
Independent. 

The  process  x f>ucn  studied  in  some  detail  since  the  mid  1960s 


but  its  history  goes  back  much  earlier;  special  cases  have  been  used  from 
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# time  to  time  since  the  early  dayi  of  stochastic  processes.  In  1915  Smolu- 

cbowski  set  up  the  following  model  for  the  number  of  particles  in  a small 
region  of  a volume  of  fluid  whose  molecules  subject  the  particles  to 
Independent  Brownian  motions.  He  assumed  that  the  number  of  particles  in 
Q has  an  equilibrium  Poisson  distribution  with  parameter  y.  The  population 
of  0 is  observed  at  time  instants  spaced  y apart  and  there  is  a probability 
P,  depending  on  the  detailed  physics  of  the  situation,  that  a given  particle  in 
Q at  ny  has  departed  by  (n-fl)T.  Independence  and  statistical  equilibrium 
then  imply  that  the  number  of  particles  diffusing  into  0 during  (nr,(n<fl)T) 
has  pgf  h(3)  s eap[-|;P(l-s)]  and  hence  if  denotes  the  particle  population 
siae  In  q at  nT,  then  ^ ^ stationary  BPI  with  f(s)  = P + (1-P)8.  Expert* 

mental  observations  have  confirmed  that  the  model  is  reliable  and  moreover 
they  gave  an  independent  estimate  of  Avogadro's  number.  The  model  was  also 
used  to  resolve  certain  paradoxes  which  arose  in  the  early  development  of 
statistical  mechanics.  We  refer  the  reader  to  [3,  ^3]  for  a detailed  account 
and  further  references.  Situations  similar  to  that  here  have  arisen  in  other 
contexts,  for  example,  in  estimating  the  speed  of  spermatoxoa  [41];  see  [20] 
[21]  for  further  discussion  and  references.  Smoluchowski's  model  is  also 
referred  to  in  an  exercise  in  [12,  Chap.  15]. 

The  BPI  has  also  been  suggested  as  a model  for  the  population  else  of 
rare  mutant  genes  or  bacteria  that  are  produced  by  an  effectively  Infinite  pop* 
ulation  of  normal  individuals  [16][50].  Another  biological  application  is  to  the 
release  of  quanta  of  transmitter  from  a nerve  terminal  under  the  stimulus  of 
an  externally  applied  electrical  pulse  [53].  This  model  is  very  similar  to 
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Smoluchowski's.  The  quanta  of  transmitter  are  viewed  as  being  contained  in 
vesicles  which  diffuse  into  the  region  of  the  nerve  terminal  and  the  numbers 
entering  the  region  between  pulse  applications  are  assumed  to  be  independent 
Poisson  variates.  When  the  field  is  applied,  each  vesicle  is  assumed  to  dls> 
charge  its  transmitter  through  the  terminal  independently  of  other  vesicles 
and  with  a constant  probability.  If  is  the  number  of  vesicles  in  the  ter- 
minal region  after  the  nth  pulse,  fX  1 is  a BPI. 

Dunne  and  Potts  [7]  proposed  the  following  model  of  a computer  con- 
trolled intersection  consisting  of  two  signalized  entries.  The  time  domain  is 
discrete  and  it  is  assumed  that  traffic  approaches  the  intersection  along  each 
entry  according  to  independent  Bernoulli  processes.  If  an  entry  has  the  green 
light,  this  remains  until  the  waiting  vehicle  queue  is  empty,  at  which  time  the 
lights  change.  Each  crossing  phase  consists  of  an  initial  fixed  period  during 
which  no  crossings  occur  and  a random  crossing  period  during  which  one  car 
crosses  the  intersection  at  each  time  unit  until  there  are  no  more  waiting 
vehicles.  The  stochastic  process  {Q^^^(n), Q^^^(n)},  where  Q^^^(n)  is  the 
number  of  waiting  vehicles  at  entry  i at  time  unit  n,  is  a two  dimensional 
random  walk.  If  is  the  time  unit  at  which  the  light  changes  to  green  for 

II 

the  nth  time  at  entry  i,  then  the  imbedded  processes  (1=1,2)  are 

BPl's.  So  also  are  where  is  the  length  of  the  crossing  portion 

of  the  nth  green  phase  at  entry  1.  The  BPI  also  occurs  as  an  imbedded  pro- 
cess in  the  one  channel  version  of  a gated  queue  arrangement  considered  by 
Leibowits  [Z6]. 

The  following  counter  model  [24]  exists  as  a functioning  device  in  the 
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Smoluchowtki's.  The  quanta  of  transmitter  are  viewed  as  being  contained  in 
vesicles  which  diffuse  into  the  region  of  the  nerve  terminal  and  the  numbers 
entering  the  region  between  pulse  applications  are  assumed  to  be  independent 
Poisson  variates.  When  the  field  is  applied,  each  vesicle  is  assumed  to  dis- 
charge its  transmitter  through  the  terminal  independently  of  other  vesicles 
and  with  a constant  probability.  If  is  the  number  of  vesicles  in  the  ter- 
minal region  after  the  nth  pulse,  fX  1 is  a BPI. 

Dunne  and  Potts  [7]  proposed  the  following  model  of  a computer  con- 
trolled intersection  consisting  of  two  signalized  entries.  The  time  domain  is 
discrete  and  it  is  assumed  that  traffic  approaches  the  intersection  along  each 
entry  according  to  independent  Bernoulli  processes.  If  an  entry  has  the  green 
light,  this  remains  until  the  waiting  vehicle  queue  is  empty,  at  which  time  the 
lights  change.  Each  crossing  phase  consists  of  an  initial  fixed  period  during 
which  no  crossings  occur  and  a random  crossing  period  during  which  one  car 
crosses  the  intersection  at  each  time  unit  until  there  are  no  more  waiting 
vehicles.  The  stochastic  process  {Q^^^(n), Q^^^n)},  where  Q^^^(n)  Is  the 
number  of  waiting  vehicles  at  entry  i at  time  unit  n,  is  a two  dimensional 

random  walk.  If  is  the  time  unit  at  which  the  light  changes  to  green  for 

n 

the  nth  time  at  entry  i,  then  the  imbedded  processes  (1^1,2)  are 

BPi's.  So  also  are  where  is  the  length  of  the  croseing  portion 

of  the  nth  green  phase  at  entry  i.  The  BPI  also  occurs  as  an  imbedded  pro- 
cess In  the  one  channel  version  of  a gated  queue  arrangement  considered  by 
Leibowits  [26]. 

The  following  counter  model  [24]  exists  as  a functioning  device  in  the 
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Electrical  Ei||sineering  Laboratoriea  at  Monasb  University.  Consider  a 
Poisson  process  and  an  independent  renewal  process,  eacb  of  whose  events 
are  registered  by  separate  counters.  The  Poisson  process  events  are  accum- 
ulated in  its  counter  in  the  normal  way  but  the  renewal  process  subtracts 
from  its  counter.  When  this  counter  empties,  the  device  emits  an  output 
pulse,  adds  r > 0 counts  to  the  accumulating  counter  and  then  interchanges 
the  input  processes  so  that  the  renewal  process  now  subtracts  counts  from 
what  was  previously  the  accumulating  counter.  If  denotes  the  total  count 
in  the  accximulating  counter  at  the  time  of  the  nth  output  pulse,  then  fX  } is 
a BPl.  In  practice  both  input  processes  are  Poisson  which  are  realized  by 
using  radioactive  sources  and  Geiger  counters.  Lee  [25]  has  given  an  exhaus- 
tive treatment  of  this  particular  case. 

Finally  we  mention  two  rather  more  theoretical  instances  of  the  BPI. 

It  has  been  suggested  in  the  dam  theory  literature  [36]  that  if  [J  1 represents 
successive  inflows  to  a dam,  then  there  are  good  reasons  for  requiring  that 
this  sequence  should  possess  the  linear  regression  property 

Most  of  the  models  used  to  realize  this  in  [36]  are  particular  BPl's.  A more 

general  analysis  was  subsequently  given  in  [30]  and  some  further  aspects  are 

considered  in  [52].  Let  fS  } be  a simple  random  walk  with  p >1/2  being  the 

probability  of  a jump  to  the  right.  We  say  that  an  over-crossing  of  height  a 

takes  place  at  time  n if  S « a,  S , « a-f  1 and  S,  * a for  some  i < n.  If 

n n-1  I , 

N(a)  is  the  total  number  of  over-crossings  of  height  a then  Dwass  [8]  has 

: . I 


s 


•hown  that  (N(0),N(1),  ■ . . ) ia  a stationary  BPl. 

The  recent  interest  in  the  BPl  began  with  Heathcote's  investigations 
[18][19]-  These  papers  seem  to  be  the  first  attempt  at  a general  treatment  in 
the  open  literature.  It  is  of  interest  to  note  that  a number  of  results  on  the 
existence  of  a limiting  distribution  and  total  progeny  of  multitype  BPX's  were 
obtained  some  twenty  years  earlier  by  Everett  and  Ulam  [11].  Our  interest 
will  center  on  the  limit  behavior  of  the  sequence  % and  we  shall  review  the 
known  results  in  the  following  section.  However,  we  shall  give  here  the  fol- 
lowing conunents  in  order  to  describe  the  new  results  in  this  paper. 

Let  m s f^(l-)  be  the  offspring  mean.  We  shall  consider  three 
cases;  m < 1 (subcritical);  m = 1 (critical);  and  1 <m  < • (supercritical). 

We  shall  defer  the  explosive  case  (m  = •)  to  a sequel.  Typically  we  try  to 
norm  the  sequence  % so  that  it  converges  in  some  sense  and  in  the  sub-  and 
supercritical  cases  it  is  known  that  this  can  be  done  If  the  following  condition 
on  If  s (hp  holds. 

Condition  A;  E(log^l)  < « 

where  lisa  random  variable  with  the  distribution  When  m < 1 the 
normlng  sequence  is  constant,  thus  is  positive  if  Condition  A holds.  If 
Condition  A fails,  then  the  norm'd  version  of  % still  converges,  but  the 
limit  distribution  is  defective.  There  is  a similar  but  slightly  more  compli- 
cated situation  for  the  critical  case;  namely  there  is  a condition  on  v ensuring 
that  S is  positive.  The  question  arises  of  discovering  something  about  the 
growth  rate  of  % when,  for  example.  Condition  A does  not  hold. 

1 
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The  author  haa  prevloualy  made  some  progress  In  this  direction 
[33][34][35].  When  m i 1 the  author  obtained  a limit  theorem  under  rather 
complicated  conditions  involving  both  v &nd  He  also 

quoted  examples  which  indicated  that  the  hypotheses  could  be  modified  to 
yield  alternative  limit  theorems.  The  main  task  of  this  paper  will  be  to  unify 
and  simplify  these  results.  We  shall  distinguish  four  different  nciodes  of 
behavior  of  the  function 

G(x)  = 1 - h(l-e‘*). 

For  all  results  that  we  obtain  we  assume  that  Condition  A does  not  hold  and 
this  is  equivalent  to  the  condition 

f G(x)dx*«.  (2) 

•*0 

The  four  modes  of  behavior  just  referred  to  are  simply  conditions  on  the  rate 


of  convergence  to  eero  of 

G(.)! 

(i) 

xG(x)— >0 

(X 

(li) 

xG(x)  — ►as  (0,») 

(X 

(iii) 

G(x)  = (x^Mx))"^ 

where  L(.)  is  slowly  varying  (sv)  at  infinity  and  either  0 < ^ < 1 or  ^*1 
and  L(x) — ►O  (x — >*); 

(iv)  as  in  (ili)  but  ^ * 0 and  L(x)— 

The  conditions  assumed  by  the  author  for  his  earlier  results 
(Theorem  4 in  [33]  and  Theorem  1 in  [35] ) are  equivalent  to  Condition  (i)  and 
the  quoted  examples  are  specific  instances  of  Conditions  (ii)  and  (iii).  The 
contribution  of  this  paper  will  be  to  provide  much  tidier  versions  of  the 
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theorems  cited  in  the  previous  references  and  with  greatly  simplified  proofs. 
We  shall  also  prove  limit  theorems  for  Conditions  (ii)  - (iv)  and  give  examples 
which,  hopefully,  will  illuminate  the  situation.  When  m = 1 we  shall  define 
a function,  also  to  be  denoted  by  G(x),  which  involves  both  v »«d  g and 
prove  limit  theorems  for  Conditions  (i)  - (iv).  The  cases  (i)  - (iii)  have  been 
dealt  with  in  [35].  The  contribution  of  the  present  paper  will  be  a simplified 
proof  for  (i)  and  a new  corollary.  Improved  versions  of  (ii)  and  (iii),  and  (iv) 
is  new.  The  present  version  of  the  limit  theorem  for  (ii).  Theorem  10  below, 
is  notable  since  it  trivially  contains  the  important  Theorem  D below. 

We  shall  now  sketch  the  nature  of  the  results  to  be  proved  and  for 

definiteness  we  shall  focus  attention  on  the  supercritical  case.  In  this  case 

Cohn  [5]  has  shown  that  if  (2)  holds  then  it  is  not  possible  to  norm  % so  as  to 

produce  a sequence  which  has  a limit  distribution  which  is  non>defective  but 

not  degenerate  at  zero.  A simpler  proof  is  presented  in  the  next  section.  The 

kernel  of  these  proofs  is  to  show  that  if  such  a norming  sequence  existed, 

denote  it  by  fc  1,  then  c /c  — ►m,  that  is  c * 0(m”^^®V  This  shows 
' n'*  n+i  n n 

that  is  growing  faster  than  some  geometric  rate.  At  this  point  it  is  worth 
pointing  out  that  another  recent  result  of  Cohn  [6]  shows  that  if  a Markov  chain 
can  be  normed  by  a sequence  fa  1 such  that  the  normed  process  has  a non* 
defective  limit  distribution,  then  necessarily  lim  sup  a ./a  <«,  that  is,  it 
can  grow  at  most  geometrically  fast.  A possible  approach  is  to  "slow  down" 
the  growth  rate  of  % by  considering  a suitable  nonlinear  function  of  and 
then  norming  the  resulting  sequence.  In  fact  the  following  version  of  this  idea 
works.  First  we  fix  ^ be  ^be  sequence  of  Seneta  constants,  to  be 
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defined  in  ^2.  In  Theorem  C below  we  point  out  that  if  (2)  holds  then 

X /c  *-  * We  shall  show  that  under  each  of  Conditions  (i)  - (iv)  there  is 
n n 

a sv  function  y^(.)  and  a sequence  such  that  { bas  a 

non-defective  limit  distribution.  In  all  cases  the  limit  distribution  is  given 
explicitly.  When  (2)  holds  we  shall  refer  to  the  "immigration  regime"  since 
the  growth  rate  of  z is  being  influenced  mainly  by  f(,  whereas  if  (2)  does  not 
hold  the  izxunigration  component  of  % serves  only  to  prevent  extinction;  the 
growth  rate  is  determined  predominantly  by  4^.  We  shall  then  refer  to  the 
"branching  regime." 

Limit  theorems  of  this  nature  occur  in  other  areas  of  probability 
theory.  The  simplest  example  occurs  when  we  take  the  logarithm  of  a product 
of  positive  independent  random  variables  and  then  invoke  the  central  limit 
theorem.  Erickson  [9]  has  proved  similar  results  for  the  forward  and  back- 
ward recurrence  timec  of  a renewal  process  whose  mean  lifetime  is  "just 


Infinite."  Finally  such  limit  theorems  occur  quite  naturally  in  the  theory  of 
explosive  Galton-Watson  process  [2]  [22]. 

The  format  of  the  paper  is  as  follows.  In  }2  we  will  review  the  knows 
results  for  the  branching  regime,  while  in  {3  we  shall  consider  the  immlgra- 
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b«en  covered  in  [27],  We  shall  then  show  how  these  results  provide  some 
Information  on  the  asymptotic  behavior  of  the  (unique)  invariant  measure. 
Finally,  in  we  shall  use  these  results  to  prove  some  limit  theorems  for 
a reversed  time  BPl.  These  complement  results  obtained  by  Esty  [10]  for 
the  non- immigration  situation. 

The  numbering  of  equations  is  sequential  within  (sub)  sections  and 
multi-place  labels  will  be  used  only  for  references  to  equations  from  other 
(sub)  sections. 

2.  The  branching  regime 

Let  [p!?^]  (n  € 19^)  denote  the  n- step  transition  matrix  of  1®* 

= Sp[j  (0  < « < i)  and  let  = f(y»))  (n«  N). 

Equation  (1.1)  implies  that 

( \ I — ^ 

P'"^s)  = (f  (s))‘  Ij  h(f  (s)).  (I) 

‘ " k=0 

If  Pq  = 0,  then  X !•  aupercritical,  'j'  • (n  — ►«»)  and  hence  all  states  are 

transient.  If  p^>0  then  [38]  S ia  irreducible,  aperiodic  and  contains 

i*  s mln[i|h^  > 0}.  It  is  of  interest  to  classify  S • When  m > 1 and  ?q>  0 


it  follows  from  (i)  that 


(n)  v.n 

p < q (h(q)) 

it* 


where  q is  the  least  positive  solutl.'in  of  £(s)  = s.  Thus  in  this  case  ^ is 
transient,  a fact  first  pointed  out  under  slightly  less  general  conditions  in  [27]. 

The  situation  is  more  complicated  when  m<  1.  Heathcote  [18][19]  proved 
that  when  m<I  and  S *1N,  then  S i®  positive  iff  Condition  A holds  and  when 
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this  i«  the  case,  the  pgf  of  the  limiting  distribution  is  given  by 


I1U)=  11  Mf.(»)). 

k=0 


(2) 


The  condition  j = IN  was  relaxed  first  by  Seneta  [43j  and  then  to  the  present 
assumption  by  Quine  [38].  We  summarize  this  as 

Theorem  A.  When  m<l,  5 ia  positive  iff  E(log^I)  < «. 

Seneta  [42]  first  pointed  out  through  an  example  that  when  m « 1 , it 
is  possible  for  S io  positive.  Foster  and  Williamson  [15]  proved  that  in 
all  cases  S i*  positive  iff 


f(s)  = s + [1-S)^^*'^{1-S) 


We  shall  be  concerned  with  the  following  specialization  of  the  critical  case: 

(4) 

where  0 < p < 1 and  «£(•}  is  sv  at  the  origin.  In  this  case  we  have  the 
following  analogue  of  Theorem  A [34]. 

Theorem  B.  Suppose  (4)  holds  and  in  addition  if  V = I then  id(s)-^» 

(s  |0)  I ds/s^{s)<«,  for  some  e > 0.  Then  S is  positive  iff 
1 

E(i*'/za  ))  < •. 

It  has  been  shown  [33]  that  when  m < 1 and  (3)  does  not  hold  then  S 
can  be  either  null  or  transient.  In  particular  if  m s 1,  * V(l-)  < • and 

y B f *(l>)/2  < m then  S null  if  o s /j/y  < 1 and  transient  if  or  >1  [28]  [29]. 

2 4 4 m 

The  case  O • I is  a boundary:  if  also  SPjJ  J <«  then  S is 

null,  and  a necessary  and  sufficient  condition  for  traussience  Is 
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m n~l 

E •*?  (-  E (l-h(£  )))  <• 

n*l  k=0 

whare  (5S][29].  These  po«slbilitie»  are  not  always  appreciated. 

For  example  In  [21]  a method  of  parameter  estimation  is  given  for  the  critical 
case  which  depends  on  the  recurrence  of  S’  It  is  erronously  asserted  that 
when  y,  h"(l-)<M,j  is  null. 

We  need  the  following  notation  to  describe  the  supercritical  ease. 

Let  r(t)  (0  <t<>logp^)  be  the  inverse  function  of  k(9)  s .logf(e  (9>0) 

and  let  r (*)  be  the  nth  functional  iterate  of  r(*)*  Fix  t.  c(0,-logq)  and 
*11  o 

let  c = 1/r  (t.);  fc  ] are  called  the  Seneta  constants.  It  is  known  [1]  that 

c c ,/c  — ►m  and  Z /c  >W  (n—>«),  where  fZ  ] is  the  sitn« 

n I nfl  n n n ^ n^ 

pie  branching  process,  Z^  = 1,  and  W is  a non-defective  random  variable 
whose  DF  has  an  atom  of  size  q at  the  origin  and  which  is  absolutely  continuouf 
on  the  set  of  positive  numbers.  Let  g(Q)  > £(0**®^).  The  following  theorem 
describes  the  growth  of 

Theorem  C.  U Xq»0,  l<m<«  then; 

(I)  X /c  V and  if  Condition  A holds,  V is  a non-defective  random 

n n ———f 

variable  with  a continuous  DF  which  is  defined  by 

E(e"®''^)  = ^ h(^(em‘“)). 
nsl 

(II)  In  addition  the  OF  of  V is  absolutely  continuous.  If  (-logh(q))/(logm) 

> 1,  V has  a bounded  continuous  density  and  if  E(l),  E(Zjlog^Z^)  < • then  V 
has  a density  which  is  continuous  on  the  set  of  positive  numbers. 
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(ill)  If  Condition  A does  not  hold  then  P(V  »«)«!.  and 

(Iv)  there  Is  no  sequence  of  constants  fe  1 such  that  a — and 

I n*  n ■—  ■ 

{X^/a^}  converges  in  law  to  any  non-defective  DF  which  is  not  degenerate 
at  the  origin. 

Parts  (i)  and  (iii)  are  due  to  Seneta  [44] [45],  (li)  to  Pakes  [32]  and 
(Iv)  to  Cohn  [5].  Seneta  [45]  also  proved  there  is  no  sequence  [a  ] such 

and  the  other  contingencies  under  (iv)  bold.  Cohn's  proof  uses 
only  the  strong  law  of  large  numbers  and  a representation  of  We  now 

give  a shorter  proof  generalizing  the  argument  used  in  [45]. 

Suppose  that  there  is  a sequence  satisfying  the  contingencies 

in  (iv)  and  let  ^0)  be  the  Laplace -Stieltjes  transform  (LST)  of  the  limiting 
DF,  F{«)-  If  * E(exp(-0X^/a^))  it  follcws  from  (1.2)  that 


Clearly 


“ IT  ))) -^^0). 

k=0  " 


Suppose  there  is  a subsequence  [n  } such  that  a , /a  — If  F (x)  is 

J n-l  n n 

A,e wwe  _ « i ^ • 


the  DF  of  X /a  then 
n n 


where  e > 0 is  a continuity  point  of  F(.). 

It  follows  from  (5)  that  ^(0)  < F(c ) < I for  all  g >0  and  where  c !• 
chosen  to  be  sufficiently  small.  This  contradicts  the  non-defectiveness  of 
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F(>).  Now  choose  a sub- seqaence,  again  denoted  by  (n^,  each  that 

a ,/a  — where  0 < (g  < m.  The  uniform  convergence  property  of  LST'$ 
nj-l  n. 

[13,  p.252]  shows  that 


j j J 

and  it  follows  that  h<f  ,(exp(-0/a  ))  converges  and  that  the  limit  is  the 

n,-i  n. 

j J 

L>ST  of  a non-defective  DF.  The  convergence  theorem  for  the  process  without 
inunlgration  and  Kinchin's  theorem  on  convergence  to  types  then  shown  that 
a mIcc  ..wkm^c  whe  re  0 < k < «».  But  we  know  that  when  Condition  A 

"j  "j-*  "j 

does  not  hold,  E(exp(-gX  /c  )) — ^0,  and  using  Kinchin's  theorem  again, 

“j  "j 

this  contradicts  the  assumed  existence  of  f* 

*•  n^ 

In  the  next  section  we  shall  prove  some  results  describing  the  growth 
of  X ^ben  the  necessary  and  sufficient  condition  in  each  of  Theorems  A-C  is 
relaxed.  When  m s 1 , one  such  result  has  been  known  for  some  timet 


Theorem  D.  U m s 1 , jj  = h'(l-),  y = f *(1-)  < m and  a « fi/y  then 

.X  . 

P(X  /ny  < X)  —>'ir(0))  r y^^*  e'^dy. 

Jq 

This  result  is  inq>ortant  since  its  conditions  will  normally  be  satisfied 
in  applications.  It  was  fouttd  independently  by  Foster  [14],  Pakea  [28]  and 
Seneta  [46] . 

3.  The  immigration  regime 

3.1  m < 1 . We  begin  our  treatment  of  the  subcritlcal  case  by  reviewing 
some  results  about  [Z^]  that  will  be  needed.  The  proofs  In  the  sequel  will 
make  repeated  use  of  the  uniform  convergence  theorem  for  regularly  varying 
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functions  and  we  shall  not  refer  to  its  use  on  each  occasion;  see  [48]  for  an 
account  of  this  and  other  properties  that  we  shall  use  about  regularly  ▼arying 
functions. 

Z 

Let  Q(s)  = lim  E(8  PS^  satisfying  the 

functional  equation  [1] 

Q(f(8))  rmQ{a)  + l-  m. 


Iteration  yields 

1 - f (s)  = «(m"(l-Q(s)))  (1) 

n 

where  w(«)  is  the  Inverse  function  of  1 - Q(l-o)  and  has  the  following  prop- 
erties [47]. 

The  function  y(x)  * x”*w(x)  is  sv  at  the  origin  and  | 1/Cy(l-)  >0  (x^O).  It 
follows  that  1 - C(l-8)  is  I -varying  at  the  origin  (equivalently,  l-Q(l-x"^)  is 
-1 -varying  at  infinity)  and  in  particular  this  gives  the  following  result  [48. 
p.l8]. 


Lemma  1.  log  (l-Q(l-s))  logbj(s)  logs  (s  ^0). 

Since  0(1 -1/x)  has  a monotone  derivative  with  reopect  to  x*^  it 
follows  [48,  pp.60,  88]  that 

(l-s)  Q'(s)/Q(s)->1  (2) 


Observing  that  <,>(•)  is  non-decreasing  it  follows  that 
h(s)/h(f^(s))<  Pj*^s)  exp(.  J logh(l-a>(m*(l-Q(s))))ds)^  1 


and  hence  if  s — ► 1 
n 


n n 
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where  ^ » 1/logm'  and  the  last  integral  arises  by  a change  of  variables  in 
the  preceding  one.  It  now  follows  from  (2)  that 


. ^ P 1 -h(»)  , 

-logp,  (v-cj  -rr«* 


■itg(i-t  (»_)) 
n n 


(3) 


= c f G(y)dy. 

*'.log(l.s^) 

n 

In  the  sequel  we  shall  always  assume  that  S positive,  that  is. 

that  (1 . 2)  holds.  The  proofs  of  the  following  theorems  all  follow  by  assuming 

one  of  Conditions  (i)  - (iv),  choosing  Ts  1 so  that  the  limit  of  the  right-hand 

' n * 

side  of  (3)  exists  and  then  invoking  the  following  proposition  [Z2l. 

Proposition  1 ■ Let  T(*)  be  a strictly  increasing  continuous  function, 

T(l)  = 0,  T(a»)  = • and  T(.)  is  sv  at  infinity.  Let  w(y)  = T'\y),  {a^)  a 
sequence  of  positive  constants  such  that  a^  — > • and  u(x)  be  a continuous 
non-decreasing  function  on  (0,c»).  If  for  a sequence  of  non-negative  random 
variables  {4^3. 


then 


E[(l-1  /w(a^x))^”3  — ► o(x) 


P(a’  T(l+g  ) < X) 


(n— X >0) 


u(x). 


Remark.  Some  changes  in  the  proof  of  Proposition  1 allow  us  to  assume  that 

T(0)  » 0 and  conclude  that  P(a  *T(e  ) < x) — ►u(x). 

n ^ n “ 

For  X ^ 1 let 


log  X 

A(x)»exp^J  C(y)dy^ 


(4) 
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end  for  0<x<l  define  y^i*}  so  that  on  [O,*)  It  is  continuous,  strictly 
increasing  and  4(0)  s 0.  E'luation  (1.2)  implies  that  4(x) — and  a change 

X j 

of  variables  shows  that  log  4(x)  « T G(logy)y  dy,  and  hence  that  4(*)  is  sv. 

‘•l 

Theorem  i . Suppose  that  m<l  and  both  (1.2)  and  Condition  (i)  hold.  Then 
[4(X^)/4(m  **))  has  the  limit  DF 

A(x)  = (1  4X)^  (x^O). 


Remark.  Condition  (i)  can  be  rewritten  as  (log  (l*a))(l-h(s)) — ^0. 

Proof.  Let  1-s  «l/A^(xA(ni  ")).  Clearly  s — ►!-  and  hence  (3)  applies. 
“ n n 

Furthermore  the  integrals  in  (3)  can  be  written  as 


= log  A((1-^„(*„))"S  - logxA(m"“) 


>logx  -f  J 


log(l-f„(8^)) 


G(y)  dy 


>logm 


Let  denote  the  integral.  Then 


where  6 


0<J  < C(-logm")[>log  (1-f  (s  ))4logm’'] 

•“  n •”  n n 

-log(l-f  (a  )) 

= 6 ( 2-^.  l)  , 

" ' -logm"  ' 

— ^0  (n — ^«*)  and  we  have  here  used  Condition  (i).  However, 


(1)  and  Lemma  show  that 


•log(i-f  (s^))  ^ .log(m"(l-Q(s^))) 
n n n 


and  'logil-Qls  ))  -log  (1-s  ) <logm’  where  in  the  last  step  we  assume 
n 11“ 

0 < X < 1 . It  follows  that  the  lim  sup  of  the  coefficient  of  0 is  finite  and 
— n 
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hence  J — >>0.  Thus  we  have  ehown  that  ) — ►A(x)  and  Proposition  1 

Q In 

•n 

(Remark)  is  now  applicable  v/ith  T(x)  = *****  *n  * *” 

It  was  shown  in  [33]  that  if  m < 1 then 

lim  n(l-h(f  ))<!(>!)  =>^  is  null  (transient).  (5) 

n 
n 

Condition  (i)  implies  that 

n(l-h(£  ))  = e(f  ) n/(-log(l.f  )) 
n n tk 

where  c(b) — (* — **>d  (*)  *****  Lemma  1 show  that  log(l-f  )<Mlogm°. 

n 

Thus  S is  null  under  the  conditions  of  Theorem  1 . 

Example  1 . The  conditions  of  Theorem  1 are  satisfied  if 

2 JL  .1 

h.  = c[j(logj)  II  log  j]  (j^J) 

^ k=2 

where  log^  j = log  j , logj^j  = *«>g(logj^  ^j),  r > 2,  J is  chosen  so  large  that 
is  well-defined  and  c is  chosen  so  that  v*  is  a probability  distribution.  It 
follows  that 

m JL  "* 

E h ~c(  II  logj)  (j->») 

lej  * ^ k=l  ^ 

and  an  Abelian  theorem  for  power  series  [13,  p.452]  yields 

_r.  I -1 

l-h(s)«<c(||  log  (1-S)“*')  (s->l-). 

'k=l  * ' 

Theorem  2.  If  m < 1 and  Condition  (li)  holds  then  {(jj/nllog  (1+X^))  has  the 
limiting  DF 

B(x)  = (x/(l  + x))*C  (x^O). 

Proof.  Let  s « 1-m^  in  (3).  Condition  (li)  Implies  that 


-log  ~ac  J 


-log(l-f 

n n 


dy/y. 


-log(l-«^) 


However,  Lemme  1 yields  log(l •£  (a  ))Mlogm  (1-Q(s  ))  end  log(l-Q(s  ))  m 

n n D n 

log(l>s  ).  Thus  If  X > 0 it  Is  now  obvious  that  ) — ►B(x)  and  Prop- 

n in 

osition  1 now  applies  with  T(x)  = (|Iog(l-fx)  and  a^  * n. 

We  remark  that  B(a|^x)  definec  an  F-dlstribution.  Condition  (ii)  and 

Lemma  1 show  that  n{l -h(f  ))  — ►ar  and  hence,  from  (5),  S is  null  if 

n ^ 

a|^  < 1 and  transient  if  a|^  > 1 . Further  details  of  the  structure  of  y are 
required  to  decide  the  case  a|^  s 1 . 

Example  2.  Let  h^  > 1 -a/log 2 >0  and 

hj  = a[(log(j+2))'^  - (log(j+3))“‘]  (jrl.Z,...). 


It  follows  that 


£ h^  a/log  j 


and  Condition  (ii)  follows  from  an  Abelian  theorem.  Conversely  it  follows 
from  a Tauberlan  theorem  for  power  series  [13,  p.4S2]  that  the  condition  of 
Theorem  2 Implies  (6).  For  this  example  S I*  null  when  a|^  > 1.  To  see 
this  observe  that  0(s)  a (1 -h(s))/(l -s)  a aS  s^/log (j-f2)  axkl  that 
a*^a(s)  - l(s)  a 0(1)  (0  < s < I)  where 


a -logs.  Moreover 

I(s)  - (T(-logT))‘*  » (t(-1o8T))"^  J(t) 
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where 


«r)=r 


(log(y>-2T))exp(.y)dy 
log(y-f2T)  - logT 


(log  (y+2T))exp(-y)dy 
(log  (y+2T))/(logr) 


and  it  can  be  shown  that  the  last  integral  converges  as  r — >*0  and  the  limit 
is  r (logy)e  ^dy  which  is  minus  one  times  Etiler's  constant.  By  patting 

•'o 

these  estimates  together  it  follows  that 

(logn)(n{l.h{f^))  - l)->0 

n 

and  hence  Bertrand's  test  shows  that  S fj”  ® *'  Thus  S i*  null. 

n^O  k=0 

Define  the  sequence  fa  1 by  G(a  ) = r/n.  If  Condition  (ill)  holds 

• tl  ^ T1  " 


-1 


then  {a  } is  ^ -varying 


-1. 


Theorem  3.  If  m < 1 and  Condition  (ill)  holds  then  fa  log(l  + X_)T  has 

1 1 ■ . I ■»  I ■ b n n ^ '■ 

the  limiting  extreme  value  DF 


C(x)  = exp(-5x'^) 


(x>0). 


Proof.  Let  s = l-exp(-xa_).  The  right-hand  term  in  (3)  can  be  written  as 


o.J, 


^log(l-f^(S^))]/[log(l.Sj] 


G(b  y)  dy 

n 


where  b « -log(l-s  ) xa  . Using  Lemma  1 it  is  not  difficult  to  show  that 

the  upper  terminal  is  of  the  form  1+  (n/|jxa^)(l+o  (1 ))  and  of  course  a/a^— ♦0. 

A mean  value  theorem  then  shows  that  -log  p/*'^(s  ) 

I n n 

and  the  theorem  readily  follows. 


i 
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Theorem  4.  U m < 1 and  Condition  (iv)  holds  then  {n[l-h(exp(-3r^))]]  has 
the  limiting  DF 

D(x)  = 1 - e”*  (x>0). 

Remark.  The  assertion  may  be  re-ejqyressed  by  stating  that  the  reciprocal 
of  the  sequence  above  has  the  limiting  extreme  value  DF  e3q>(>x~^). 

Proof.  Let  l-s  = [l+ 1 /(.logh'^l-(l-h-)/(l+xn(UhJ))]"* 
n u o 

^ -logh‘Nl-(I-ho)/(l+xn(l-hQ))). 

Since  G(.i  is  sv  at  infinity,  x®C(x)— (x — ►«»)  for  any  6>0  and  this 
allows  us  to  infer  that  >log  (1-s  ) > K n^  , 0 < K < •,  A*  > 0 and  hence  that 

[log(l.ysJ)]/[log(l.s^)]  = 1 - (n/(;log(l.s^mi+o{l)) 

and  n/log  (l>s  ) — ^0.  By  proceeding  in  a similar  manner  to  the  proof  of 
11 

Theorem  3,  we  obtain 

.logp/"\s^)  ^nG(.log(1.8^)) 
in  n 

~nG(-log(l.h"V-(l-hQ)/(l+xn(l.hQ))))] 

S n(UhQ)/(l+xn(l-hQ))--^x‘*. 

Proposition  1 (Remark)  may  now  be  applied  with  T(x)  « (l»h(e*^*))”^-  (l-h^l’^ 

a « n and  u(x)  > 9xp(-x  xod  the  theorem  follows, 
n 

Under  the  conditions  of  Theorems  3 and  4, 

l.h(fj  = G(-log(l.f^))  ~(n/u)"VL(n) 
where  we  have  used  (1 ) and  Lemma  1 and  it  follows  from  (5)  that  S 


transient. 
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Example  3.  Let 


= cj'^log  j)  * 


(j=2,  3,  . . . ; 1 < a < 2) 


where  c is  a constant.  It  was  shown  in  [33]  that 

1 - h(8)  ~ c(a-l)(-log(l-s))^’® 

and  it  follows  that  the  conditions  of  Theorem  3 are  fulfilled  with  a = 

( -1)"^ 

Moreover,  it  is  clear  that  we  can  replace  a^  by  a^  = c(a-l)(n/u)  ® since 

a ~ a'  . 
n n 

Example  4.  '’^e  modify  the  previous  example  by  letting 


k*  1 


for  j>J,  *■^2  and  b > 0.  It  follows  that 


and  hence 


Thus 


1-h.  ~ (O'b)(log^j)' 


l-h(s)  ~ (c/b)(log^(l-s)’S"^ 


G(x)  ~ (c/b)(log^_^x)' 


(8-+-U). 


and  this  satisfies  Condition  Iv.  Furthermore  since  X we  may  replace 

n 

l>h(exp(-l  /X  ))  by  (c/bKlog"*^  X )'^  in  the  assertion  of  Theorem  4 and  the 
n r n 

remark  then  shows  that  fn‘*  ^^(log***  X )1  has  the  limiting  extreme  value  DF. 

' r n 

exp[-(c/b)x  (x  > 0). 

It  is  interesting  to  note  the  role  of  Lemma  1 in  each  of  the  proofs 
above,  namely,  that  because  of  it  the  offspring  distribution  figures  in  the 
results  only  through  the  offspring  mean  m.  We  can  make  some  comparison 
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of  the  growth  rate  of  % for  each  of  Theorems  1-4.  First  note  that  the  tail  of 

ff  grows  longer  as  we  go  from  Conditions  (i)-(iv)  and  hence  we  expect  the 

growth  rate  to  increase  as  we  go  from  (i)-(lv).  This  is  clearly  so  for  Theorems 

2 and  3 since  after  normalising  by  log(.)  a linear  norming  v/orks  in  Theorem 

2,  but  the  norming  sequence  fa  1 in  Theorem  3 grows  at  moat  algebraically 

fast,  but  a /n — If  Condition  (iv)  is  in  force  then  l-h(expx  S = 
n 

G(-log  {1 -exp  x"^ ))  ~ G(log  (1 + x )).  Thus  we  can  always  re-express  Theorem  4 
as:  {n’^L  (log  (1+X^))}  has  the  limit  DF  exp(-x  ^).  In  this  case  the  non- 
linear norming  Is  much  more  severe  than  for  Theorems  2 and  3.  We  have  the 
same  situation  for  Theorem  1 since  under  Condition  (i)  it  is  easily  shown  that 
y^(e  ) is  sv  at  infinity.  However,  it  is  also  true  that  the  norming  sequence 
(A(m  ^)}  Is  sv.  It  does  not  seem  easy  to  exhibit  specific  examples  of 
However,  we  can  derive  bounds  in  the  case  of  Example  1 as  follows.  It  fol- 
lows from  the  definition  that  xj((x)/\(x)  ^ Kh(l-x  S and  hence  for  Example 
1 we  obtain 

r -1 

xA'(3c)/A(x)  ^ c ^ ■ 

Let  c > 0 be  given.  If  x^  is  large  enough  it  follows  that 

c(l-e)(log^^jX  - log^^jX^)  <log(A(x)/A(xQ)) 

<c(l+c)(log^^jX  - 

or 

Kj(log^x)*^^^"  < a(x)  < K2(log^x)‘^^*'*’®^  (7) 

Cohn  [4]  has  given  a theorem  which  gives  conditions  which  permit 
almost  sure  convergence  of  an  in- homogeneous  Markov  chain  to  be  inferred 

( 
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from  its  convergence  in  law.  He  showed  that  this  was  applicable  to  the  chain 

fT(l+Z  )/a  I where  fZ  1 is  an  explosive  branchinj;  process,  T(.)  is  a 
^ n n ' *•  n 

certain  sv  function  and  * certain  sequence  of  constants.  Despite  the 

similarity  of  this  situation  and  those  considered  in  this  paper,  convergence  in 
law  is  the  most  we  can  assert  for  our  results.  This  follows  because  each  of 
the  sequences  in  Theorems  1-4  is  mixing  [40,  Theorem  2j  and  hence  cannot 
even  converge  in  probability  [39]. 

3.2.  1 < m < »■  Recall  the  notation  introduced  before  Theorem  C and  let 

K(0)  = -log  ).  The  following  facts  are  demonstrated  in  [47]; 

K(0)  = 0X(9) 

where  \\»)  Is  strictly  decreasing,  sv  at  the  origin  and  ^EW  (9  ^0).  Let 
u(»)  be  the  inverse  function  of  K(.).  Then 

T (t)  = K(m""w{t))  (0<t<-logq)  (I) 

n 

and  it  is  clear  that 

K(l)  = tjj,  ult^)  = 1 (2) 

and  hence  that 

u(p^)=m'“,  K(m'”)  = (3) 

where  inverse  function  of  0(»);  it  is  defined  for 

q < 8 < 1 . Clearly 

K(0)  = -log0(0),  w(t)  = ^(e'*)  (4) 

and  hence  the  regular  variation  properties  of  K(«)  and  &>(•)  are  inherited 
by  1 -0(0)  and  ^»)(l-8),  respectively.  In  particular 

(l-a)^'(8)/^(s)->-l  (a-*-!-). 


(5) 
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i* 

Finally,  alnce  ^(0  satisfies  th«  functional  equation  ^(Om)  3f(0(0))  [1],  it 
follows  that 

f^(a)  » 0{m“0(s))  (q<s<l).  (6) 

The  following  lemma  Is  a consequence  of  the  regular  variation  property  [48, 
p.l8j. 

Lenama  2.  logK(O)  ~ log(l  - 0(0))  ~ log0  (0  ^ 0) 

logw(t)  logt  (t  ^0),  log0(s)  log(l>s)  (s  — >1-  ). 

The  results  on  the  convergence  of  menMoned  In  ^ tell  us 

that  f (s®»)  — ►0(-log  s)  (0<s<l)  and  hence  If  a — £ (s^®)“^l. 

Thus  it  follows  from  (2.1)  that  p/"Ns^)  *,pi“^(s^).  If  I (s)  » 

i n o n n 

r”  X 

exp  I logh(0(m  0(s)))dx  It  follows  using  (6)  that 

1 ^ p‘“\s)/I  (s)  s h(s)/h(f  (a)) 

^ o n n 

and  hence  if  s -►I-,  P,J”^  (s  I (s^").  Changing  variables  in  X (s) 

n o n n n a 

and  using  (5)  we  eventually  obtain 

-log(l-s^n) 

-log  (s^^“)  ^ C J C{y)dy  (7) 

-logd-f  Js^^n)) 

n n 

where  > 1/logm. 

We  are  now  ready  to  obtain  theorems  for  > Pn^n  analogues  of 

Theorems  1-4.  For  the  next  theorem  we  define  exactly  as  in  (3.1.  It 

j was  nxentioned  in  (3. 1 that  Condition  (i)  implies  that  infinity. 

I We  shall  need  to  use  this  fact  in  the  proof  of  the  following  theorem. 
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Theorem  S.  U (i.2)  and  Condition  (i)  hold,  the  eequence 

{4(Y^)//^(m“))  ha»  the  limiting  DF  A(x). 

Proof.  Let  (1-s^)"^  = A \xA(m°))-  The  integral  occurring  in  (7)  can  be 
written  as 

-log(l-8^") 

r G{y)dy  - log[y^((l-£^(s^  ))  “*  (8) 

logm 

By  using  the  mean  value  theorem  we  can  write 

Pn 

-log  (1-s  ) s -log  6 -logP  - log(l-s  ) (9) 

n Q o n 

where  0 — ^1  (n — It  follows  from  (3)  and  Lemma  2 that 
n 

(-logp  )/(logm") —►1.  (10) 

n 

If  0 < X < 1 , -log  (1  -s  ) < logm*^  and  hence  it  follows  as  in  the  proof  of 
••  n •• 

Theorem  1 that  the  integral  at  (8)  tends  to  zero.  To  proceed  further  we  prove 
Lemma  3.  Assuming  only  that  s^ — ►!-,  logm  ) »»*log{l-»jj) 

Proof.  First  observe  that  (4)  and  the  exponentiated  form  of  (10)  show  that 

P P 

a(8  ")«mw(1-s  )«»w(p(l-s  ))  and  hence  from  (3) 

~ a n n n 

m"A(s^")  (4p  (1-8  ))/u(p  ) » (!-■-)  y(P_(l-*))/y(^_) 

^ n n o n ii'O  ii*ii 

where  y(t)  * t'^b)(t)  is  sv  at  the  origin.  Since  w(x~S  has  a monotone  deriv- 

ative with  respect  to  x**  it  follows  that  [48,  pp.  60,88] 
where  e(t)— ^0  (t^O).  Integration  yields  the  representation 

1 *0  -1 
y(t)  * t“  exp  ( - J e(y)y‘  dy) 
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whence 

a . 

(*  | «(y)y’  ^y) 

n^ 

= exp  (fi„log(l-8_)) 

'A  D 

where  fl  — >*0  and  we  have  used  a mean  value  theorem, 
n 

It  now  follows  that  [logm"i5(s^  )]/[log  (1 -s  )]  ^ 1 4-  t?  and  the  assertion 

n n n 

follows . 

Returning  to  the  proof  of  Theorem  5,  observe  that  Lemma  3 implies 
n P X 

that  m — ^0  and  hence,  if  R(x)  = A(e  ) it  follows  from  (6)  and  Lemmas 

2 uid  3 that 

\ R.[-log(l-f  (s^))]  ~ R[-logm"0(s%] 

~ R(-log(l-s^)J  = = *A(*"°)- 

(n)  Pb  ^ 

The  proof  is  now  easily  completed  since  Pj  '"jj  ) ~ 

the  slow  variation  of  R(x)  was  used  in  the  last  sequence  of  steps. 

Theorem  6.  U 1 <m<ai  and  Condition  {iV  holds,  then  [(jj/n/ log(lf 
has  the  limiting  DF  (B(x). 

Proof.  If  1 - s = m”***  (x  > 0)  it  follows  that 

— n 

p 

Y log(l-s  ") 

>logE(s  ")#«arlogf 2 * 

**  ' Pn  ' 

log(l.f  (s  "» 
n n 

Y 

It  follows  from  (6),  Lemmas  2 and  3,  (9)  and  (10)  that  E(s  ** ) — ^B(x),  whence  | 

n f 

I the  assertion.  I 

i ! 
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Let  Condition  (ill)  hold,  define  a by  G(a  ) = r/n,  let 

n n ^ 


P * n 

= -log(I-8^*')  and  = « • If  follow*  from  (10)  that  (logjj^a^ — ►O 

and  hence  from  (9),  b xa  . Thus  we  finally  obtain  G(b  ) ~ x ^(r/n).  We 

n n n * 


can  rewrite  (7)  as 


.logE(8^")  --(J  J 


[log(l-f  (s^"))]/[log(l-s^")] 


C(b  y)dy. 
11 


Denote  the  lower  terminal  by  t^.  Then  since  1 - ^(g)  = 9)i(0)  and  ^(s)  = 
(f-*)g(l-B)  where  )x(«)  and  g(>)  are  sv  at  the  origin,  it  follows  from  (6) 


P P P 

log(l-f  (8  “))  s logm”  + log(l-s  “)  + logg(l-s  ") 


n ^n 

+ logn(m  ^(s  )). 


This  yields 


" Clog(l-8^n) 


Pn 

Clogg(l-8^  ) 
nlog  (l-B^^n ) 


n P 


n , 


logji(m  ^(s^  ))  (logm  ,[,(8^  ) 


n , .'^n. 


logm  a(8  ) 

^ n 


1 /I 

nlog(l-s  ) 
n 


Clearly  the  second  term  in  parenthesis  — >-0  [48,  p>18]  and  also  using 

Lenama  3,  it  follows  that  the  third  term  also  — >0.  Finally  using  (9)  and  (10) 

we  obtain  1-t  « (n/rxa  )(l-fO(l)).  We  can  now  proceed  exactly  as  in  the 
n ^ n 

proof  of  Theorem  3 to  obtain 

Theorem  1.  U 1 <m  < m and  Condition  (iii)  holds  then  fa‘^log(l-f  Y ))  has 


I 


i 

I 

\ 
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the  limiting  DF  C(x). 

Theorem  8.  If  1 < m < » and  Condition  (iv)  holds  then  fn(l-h(cxp  - Y”^))l 
has  »he  limiting  DF  D(x). 

Pn 

Proof.  The  relation  (11)  still  holds  if  b = -log(l-s  ) and 

n n 

1-8^  = [l+l/[-logh’^l  - (I.hQ)/(l  + xn(l-hQ)))ir*. 

Using  the  last  part  of  the  proof  of  Theorem  7,  we  can  proceed  similarly  to  the 

proof  of  Theorem  4 once  we  show  that  b ~ log(l-8  ).  To  do  this  it  suffices 

n !i 

to  show  that  (-log(l-s  ))/n — 

n 

Let  H(.)  be  the  inverse  function  of  C(*);  It  Is  defined  on  (O.l-hg]. 

Since  H(y)  = -log  (1 -h'^1 -y))  we  have  that  -log(I-s  ) ~H((l-h^)/(l+xn(l-h^))) 

no  0 

and  hence  it  suffices  to  show  that  n*^H(n”S“^«.  Let  0<e<l  be  given. 

Then  for  any  n large  enough,  n > N(e)  say,  then  l-e<  C(e~^n}/C(n)  and 

(I-e)  <nG(tt).  These  imply  that  1 < nG(e  ^n)  and  hence  that  n”^H(n”^)  > 

The  form  of  Theorems  5-8  makes  it  clear  that  the  examples  and 

discussion  in  {3  1 apply  here  with  obvious  changes.  Finally  we  mention  that 

if  E(Z,  log^Z.)  <•  then  we  can  replace  Y by  m”"x  in  Theorems  5-8.  • 

* t n n 

3.  3.  m g 1 . 

Let  U(s)  (0  < s < 1)  generate  tl>e  invariant  measure  of  the  process 
(Z^)»  normalized  so  that  U(Pq)  = 1.  Suppose  now  that  (2.4)  holds.  It  is 
known  [Sl]  that 

U(s)  = (,;(1.S)»'M(I-S)]‘‘  (1) 

where  M(s)m.|£(s)  (s^O).  The  function  U(l-s)  has  an  inverse  g(»)  defined 
on  (0,»)  and  hence 
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- \ Iv 

g(x)  = X K(x)  (2) 

where  K(«)  is  sv  at  infinity.  The  functional  equation  U(f(s))B  HU(8) 

[1 , p.  68]  may  be  iterated  to  yield 


1 - f («)  = g(n  + U(«)) 


and  this  may  be  used  [34]  to  show  that  if  s — > i- 

T\ 

/ , U(f  (s  )) 

-logp/’  is  r " “ G(y)  dy 


=I 


n+  U(s  ) 
n 


U{s  ) 
n 


(3) 


(4) 


where  we  now  define 


G(x)  = 1 - h{l-g(x)). 

These  relations  reflect  the  fact  that  the  results  to  be  derived  depend  on  a more 

intimate  interaction  between  V 9 than  is  the  case  for  m ^ 1 . The  roles 

played  by  log(«)  and  e * when  m 1 are  taken  here  by  U{.)  and  g(»). 

When  m < 1 we  made  frequent  use  of  the  relation  -logCl-f  (s  ))  ~a/r  - 

n n * 

log(l>Q(8  ))  and  a similar  one  was  used  when  m > 1.  In  the  critical  ease  the 
n 

analogous  relation  is  the  identity  U(f  (s  ))  = n 4 U(s  ).  Thus  there  are  manr 

n n n 

similarities  in  the  structure  of  the  critical  case  and  that  of  the  non>eritlcal 

cases.  However,  there  are  also  great  differences  as  our  results  will  show. 

We  shall  assume  in  the  sequel  that  ^ is  not  positive,  that  is 

f G(y)dy  = •. 

’'O 

Let 


(5) 
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U(l-I/x) 

y^(x)  = exp  I G(y)dy  (x^fl) 

and  extend  the  definition  to  [O,'*)  so  that  /.(.)  is  continuous,  strictly 
Increasing  and  ~ 

Theorem  9.  l.,et  m = 1 and  (5)  and  Condition  ^1)  b^ld.  lf^a^=(l-f^)  ^ then 
has  a vnirorm  digtrihu' on  [0,1]. 

Prr.cf.  Let  (!-3  )*^  - a'^  (-:*;(■'  It  follows  as  in  the  proof  of  Theorem  1 , 

^ A ***  P 

that 

, , .n+  U(s  ) 

-log Fj“'  (s^) log  X + J G(y) dy 

*^n 

and  since  if  x < 1 U(s  ) < U(f  ) = n,  the  integral  is  dominated  by  nG(n)  — >0. 

n ~ n 

The  assertion  readily  follows. 

There  is  a difference  between  this  result  and  the  corresponding 
resxilts  for  m 1 , namely  that  /^l.)  is  sv  at  infinity  only  as  a consequence  of 
the  hypotheses  of  the  theorem.  To  sec  this  observe  that 

xA'W/aW  = [x'V(l-l/x)/U(l-l/x)]  U(l-l/x)G(U(l.l/x)I 


and  the  term  in  square  brackets  — >%/  (x — It  is  not  true,  without  assump- 
tions, that  ^ were  the  case  then  the  norming  sequence  in 

Theorem  9 could  be  replaced  by  I*  argue  directly  as 

foUows.  Suppose  we  let  = A"^  (x  above.  We 

would  then  be  led  to  the  Integral 


n+U(l-n‘^^*') 


-1/v 

C(y)dy  < nG(U(l-n  ' )) 
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and  since,  by  (1),  n/U(l = (VM  ‘ the  following  coroUary  ia 

apparent. 

Corollary.  If  the  conditions  of  Theorem  9 hold  and  if  lim  inf  /(s)  >0  then 
(/^(X  has  a limiting  uniforn^  distribution  on  (0,1). 

If  s’®(l-h(l-8))  is  8V  at  the  origin  it  ia  readily  checked  that  the 
conditions  of  Theorem  9 can  hold  only  if  6 - V.  If 

for  j sufficiently  large  and  r > 1 then  if  0 < t/  < 1 


l.h(l-e)  (c/v)r(»-«^)8’^(  ) 


• l/V 

Using  (2),  logg(x)~logx  (x--><»)  and  hence 


— 1/  ' 

:G{x)  ^ (c/l/)r(l-»')(K(x))*'  ( II  log  X 

'‘k=l  ' 


Thus  the  conditions  of  Theorem  9 will  certainly  be  fulfilled  if  K(x)  is  bounded 

1-f  V 

and  this  is  the  case  when,  for  example,  f(s)  = s + (b/(l+ i;))(l-s) 

(0  < b < 1,  1 < P < 2).  This  example  satisfies  inequality  (3.1.  7)  provided  c is 
replaced  by  (cP/b)(l+ P)  r(l -P).  Moreover  the  conditions  of  the  corollary  are 
satisfied,  but  It  is  evident  that  if  r = I,  A(e*)  is  not  sv. 

A further  example  was  constructed  in  [35]  for  which  ^ = h'(l-)  < • and 
f(s)  - s ~ K(1-8)^  (log  (l-s)"S^”®  where  0<K<*and0<c<l.  It  can  be 
shown  that  for  this  example 

expI((l-t)^/Kc)(logx)®J  < A(x)  < *xp[((l+  e)^/Kc)(logx)®  ] 
so  that  /^(,)  increases  much  more  rapidly  than  for  the  previous  example.  If 
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V = 1 the  extra  condition  in  the  corollary  ia  alwaya  aatiafied  aince 

/(1-a)  = (l-a)'*  (1  -£(a)) 

where  f (a)  s (i.f(8))/(l.a)  ia  a pgf  and  hence  i£(I-a)  geaerateathetailprob- 
abUitiea  of  the  diatribtttion  defined  by  T(.).  Thua  ^(x)  increaaea  aa  x^O. 

On  the  other  hand,  it  ia  not  difficult  to  construct  examplea  for  which  the  extra 
condition  ia  not  aatiafied. 

Although  the  two  following  theorema  luive  a quite  different  character 
from  our  other  reaulta,  they  are  in  a aenae,  dual  to  Theorema  2 and  6,  and  3 
«nd  7,  reapectively. 

Theorem  10.  U m > 1 and  Condition  (ii>  holda  then  [(l-f^)X  ) haa  a 
a limiting  diatribution  whoae  LSTje  (1  + g*^)”* 

Proof.  Let  ■ exp(-9  (1-f^))  and  uae  (4)  and  note  that  (1)  implies  that 

U(a  J ^ e‘*'n. 
n 

By  ualng  (2)  we  aee  that  the  hypotheaea  are  equivalent  to  (2.4)  and 

l.h(a)  .,av(l-a)*'^(l-a)  (a-wl-). 

Kawasu  and  Watanabe  (23]  started  with  these  assumptions  obtained  the 
eoncluaion  of  Theorem  10  by  ualng  a much  more  complicated  metted.  Finally 
if  y<«  then  V*l,  l-f^«.{ay)  * [1]  and  since  i£(e)— (a  — >0)  It  follows 
*****  f(*)  ^ (*y)  • lf»  in  addition,  fi  < ^ it  follows  that  Condition  (ii)  la  sat- 

sfied  with  a > o and  hence  we  obtain  Theorem  D.  This  waa  shown  la  a 
nore  roundabout  way  ia  (34).  It  waa  also  shown  la  this  reference  that  S 
• null  when  a < 1 and  transient  when  a > 1. 
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Theorem  1 1 . Let  m « I,  Condition  (lil)  hold  and  let  »^  = I-h  '(1-1/n).  Then 
hne  the  limiting  sUble  law  whose  LST  ie 

Proof.  Let  B^*exp(-0aj  and  Now 

n » 

• ince  U(h‘*(l-x))  * g”^(x).  It  follows  frcm  Condition  (iii)  that  ^(x)  is 

sv  at  the  origin  and  hence  that  n/b^  ->0.  Condition  (III)  and  (7)  also  show  that 

G(b  ) ~ 0^‘'n“*.  Equation  (4)  can  bo  rewritten  as 

1 + n/b 

-log  (s  ) ~ b f G(b  y)  dy  ~ nG(b  ) ->0^*^ 

I n njj  n n 

and  the  theorem  follows. 

The  hypotheses  of  the  theorem  are  equivalent  to  the  requlrementsr 
Eq.  (2.  4)  and 

l-h(s)  = (l-o)*p(l-s)  (0<6<l')  (8) 

where  |i(.)  Is  sv  at  the  origin  and  if  8*2,  p(s)/|£(s) — ►*  (a^O).  For  it  Is 
easily  seen  that  Condition  (ill)  and  (2)  imply  (8)  with  6*  and  |i(s)  ^ 
(p/{o))^/L((l-s)"*'/|^(s)),  and  conversely  (2)  and  (8)  imply  C ndltlon  (ill)  with 
t*t/v.  When  y<*  Theorem  11  was  given  in  [33,  Theorem  8,  Remark  2].  We 
see  that  the  Conditions  of  Theorems  10  and  11  require  |f  to  have  a regularly 
varying  tail  with  index  in  0)  and  this  oosaibility  is  also  allowed  by  Theorem 
9.  In  particular  the  hypotheses  of  Theorems  9’11  imply  that  possesses 
moments  of  some  algebraic  order  < 1 . A boundary  case,  not  yet  covered, 
occurs  if  possesses  no  algebraic  moments.  This  occurs  If  has  a sv  tall 
and  the  following  theorem  covers  this  contingency. 


i 


I 
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Theorem  12.  U m = 1 and  Condition  (iv)  hold  a then  {n(l-h(exp(-X^^)))}  has 
the  limiting  DF  D(x). 

Remark.  The  hypotheses  Imply  that  (8)  holds  with  6*0  and  conversely.  If 
(Z)  and  (8),  with  6*0.  hold  then  so  does  Condition  (iv). 

Proof.  Let  s = h’^  (I  - (I -h„)/(l+ xn{l -h^)))  and  b = U{8  ).  Using  the 
' ■ n 0 0 n n 

argument  at  the  end  of  the  proof  of  Theorem  8 it  follows  that  n/h^  —>0. 

Furthermore,  G(b  ) = h(s  ) ~ (xn)"^  and  hence,  by  proceeding  exactly  as  in 
n n 

the  last  few  lines  of  the  proof  of  Theorem  11,  we  find  that  — ►exp(-x’^). 

The  assertion  follows  from  Proposition  1. 

Equation  (3)  implies  that  l-h(f^)  = G(n)  and  hence  it  follows  from 
(3.1.  5)  that  S transient  under  the  conditions  of  Theorems  11  and  12.  The 
remark  shows  that  all  the  examples  of  v in  $3. 1 satisfy  the  hypotheses  of 
Theorem  12.  Thus  it  follows  tluit  for  Examples  1>4,  the  following  sequences 
converge  in  law  to  the  extreme  value  OP  exp(-x  ^): 

-1  — 

Example  1 (r>2):  {(cn)  ||  log  X } 

k=l 

Example  2t  {(an)  ^ logX^] 

Example  3:  ((cn(6-l))  * (logX^)^  (0<6-l  < 1) 

Example  4:  {(b/cn)  b >0). 

Many  other  examples  could  be  constructed. 

4.  The  Markov  branching  process  allowing  immigration. 

In  this  section  wc  shall  consider  the  Markov  branching  process 

I 

i 
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allowing  immigration,  denoted  by  [X^;  t > 0}.  The  intuitive  picture  for  this 
process  is  as  for  the  BPI  except  that  each  individual  has  a lifetime  that  has  an 
exponential  distribution  with  mean  \ , the  lnter>immigration  times  are 

exponentially  distributed  with  mean  6 ^ and  these  times  are  mutually  inde> 
pendent  and  independent  of  the  other  elements  of  the  process.  The  number  of 
individuals  living  at  time  t is  X^.  Let  a(s)  = b(s) » 6(h(s)-l ) and 

F(8,t)  is  the  pgf  of  the  process  without  immigration  descended  from  a single 
individual.  It  is  known  that  [l] 

t = dy/a(y). 

•*8 

We  assume  that  F(I-,t)  « 1,  that  is,  dy/a(y)  = • for  each  sufficiently 

•’l-e 

small  e • 

It  is  known  from  [49]  that 

X t 

Pj(s,t)  = E(8  *|Xg=  i)  = (F{8,t))  exp(J  b(F(8,T))dT)  (1) 

A change  of  variables  yields  Pjj(8,t)  = ^[{s)/  ||  (F(s,t))  where  1^(3)  = 

8 

exp(>  I (b(y)/a(y))dy)  is  the  generating  function  of  the  stationary  measure  of 

[X^}  [31]  [54].  The  stationary  measure  is  unique  up  to  constant  multipliers. 
Since  Theorems  A and  B are  consequences  of  (2.3),  it  immediately  follows 
from  the  form  of  [{(O  that  with  a trivial  rewording.  Theorems  A and  B 
apply  to  {X^). 

We  shall  work  directly  with  (1)  and  deduce  asymptotic  expressions 
similar  to  (3.1.3),  (3.2.7)  and  (3.3.4).  Once  this  is  done  the  proofs  used  in 
)3  carry  through  here  with  only  trite  changes.  This  approach  seems  to  be 
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more  useful  here  than  the  embedding  method  that  has  been  used  for  limit 
theorems  describing  the  branching  regime  [31],  Let  a = a'(l-)  =x(ni-l)  < 


g < 0.  Let  [Z^)  denote  the  Markov  branching  process  and  =1.  It  is 


known  that  E(s  ^|Z^>0) — ►Q(s)  and  that 

1 -s 
'0 


A 9 

1 - Q(l-s)  = exp^cr  J dy/a(y)^ 


1-8 

8 exp  J (14 


•t] 


xa(l-x  S 


-1. 


Since  (l  + a/xa{l-x  )) — ►O  (x — ►*)  it  is  clear  that  the  coefficient  of  s Is  sv 
at  the  origin  and  hence  both  1 - Q(l-8)  and  its  inverse  u(«)  are  1-varying  at 
the  origin.  The  function  Q(<)  satisfies  the  functional  equation  l-Q(F(s.t))  = 
e^(l-Q(s))  whence 

1 - F(s,t)  = «(e"‘(I-Q(s)) 
and  it  follows  that  If  s^  — ► 1 - then 


.log(l-F(8^.t)) 


-logP  (8  ,t)  ~ 5 r C(y)dy 

•’-log(l-s^) 


where  = 6/\(l-m)  and  C(x)  = l-h(l-e  ).  Lemma  1 continues  to  apply  and 
we  now  see  that  we  have  virtually  the  same  situation  as  in  §3. 1.  Define  /^(O 
as  at  (3.1.4).  The  following  result  is  obvious. 

Theorem  13.  Let  g < 0,  I*  G(y)dy  = <■  and  = fi/xCl-in)-  If  Condition  (i) 

•*0 

holds  then  ^ ))  has  the  limiting  DF  A(x). 

If  Condition  (ii)  holds  then  ((it  (l-m))*^log  (14XJ1  has  the  limiting  DF  B(x). 
If  Condition  (iii)  holds  and  a(t)  is  defined  by  G(a(t))  * -t/g  then 


I 
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[(1 /»(t))  log  (1 +Xj))  haa  the  limiting  DF  C(x). 

If  Condition  (iv)  holds  t><gn  [ - tb(exp  (-X^*))]  has  the  limiting  DF  D(x). 

. a 

0 < g < a».  Let  k(0,t)  = -logF(e  , t)  and  denote  its  inverse  function  by 

T (x,  t):  it  exists  if  0 < x < -log  F(0,  t).  Choose  XQe(0,-logq)  and  let  c(t)  = 

l/^fx^.t).  It  is  known  [3lJ  and  references  therein,  that  c{t) 

c{t+ T)/c(t) — and  Z^/c{t)^^^^  W where  W is  a non-defective  random 

variable  whose  OF  has  an  atom  of  size  q at  the  origin  and  is  absolutely  con- 

-flW 

tinucus  on  the  set  of  positive  numbers.  Let  ^{0)  = E(e  ” ).  Then 

0(0)  = F(0(0e"^^,t) 

and  if  ^(s)  is  the  inverse  function  of  ^(.)  it  is  shown  in  [31]  that 

s 

0(s)  = A(l-8)exp  ^ J (€<•<)) 

e 

where  A = 0(e)/(l-c)  If  follows  that  ^(1-s)  is  1 -varying  at  the  origin  and 
so  also  is  (1-0(0).  Defining  K(.)  and  u(«)  as  in  §3.2  we  see  that  Lemma  2 
holds. 


Returning  to  the  immigration  process,  it  was  shown  in  [31]  that 
X^/c(t)^^  V where,  if  Condition  A is  satisfied,  V is  a non-defective 
random  variable  with  an  absolutely  continuous  distribution  and 


-oV  “^0 

E(e  ) = (0(0))  exp  I (b(y)/a(y))dy. 

•'l 

If  Condition  A is  not  satisfied  V •,  and  we  now  assume  this  to  be  the  case. 
Lot  p(t)  * l/c(t),  = p(t)X^  and  s^ — ►!-.  It  is  easily  shown  that 

Y -logd-i^^ 

-logE(s  r G{y)dy 

’'-log(l-F(sJ:t)) 
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where  now  = fl/x(m-l).  It  should  now  be  clear  that  we  can  write  down  a 
list  of  results  with  almost  the  same  wording  as  Theorem  13;  replace  "a  < 0" 
by  ”0  < o < •, " "X^"  by  " Y^,  " and  in  the  statement  corresponding  to 
Condition  (i)  replace  "A(® 

g g 0.  It  is  known  [17]  that 

s j 

U(s)  = r (f(x)  -x)  dx 
•*0 

generates  a stationary  measure  for  {Z^}  ‘bat  U(F(8,t))  *X*+^t*) 
hence  that 

1-F(8,t)  = g(xt  + U(8)) 

where  g(«)  is  again  the  inverse  of  Assuming  that  (2.4)  holds  here, 

it  immediately  follows  that  (3.  3. 1)  and  (3.  3.2)  apply  in  the  present  context. 

Again  defining  G(x)  = 1 - h(l-g(x)),  (1)  can  be  transformed  to 

Xt+U(s) 

-logP  (s.t)  g (6/x)  G(y)dy. 

•'u(s) 

The  proofs  of  (3. 3 apply  here  but  there  are  some  changes  in  parametrlzations. 

m 

Theorem  14.  Let  g = 0,  I*  G{y)dy  = •»  and  (2.4)  hold. 

•'o 

If  Condition  (i)  holds  and  a(t)  g (UF(O.t))*^  then  {a(Xj)/A(«(*)>)  bas  the 

limiting  DF  (1  If.  in  addition,  lim  inf  i£(s)>0,  we  can  replace 

1/1/ 

A(*(t))  with  A(t  )• 

If  Condition  (ii)  holds  then  [(l-F(0,t))X^)  has  the  limiting  DF  whose  LST  is 

If  Condition  (iii)  holds  and  a(t)  = I - h‘^(l-t"S  (*^1)  then  {a(t)X^}  has  the 
limiting  stable  law  whose  LST  is  exp  (-68*^^). 

If  Condition  (iv)  holds  then  {t(l-h(exp(-X*S)J)  has  the  limiting  DF  D(6x). 
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fi,y  it  was  shown  in  [31]  that  [X^/t]  has  a limltinf  ganima 
law.  This  follows  ixnmediately  from  the  Condition  (ii)  statement  above. 

5.  Asymptotic  behavior  of  transition  probabilities  and  the  stationary  measure. 


A not  inslgnificcnt  portion  of  the  literature  on  the  simple  branching 
process  has  been  devoted  to  the  determination  of  the  rate  of  decay  of  its  n>step 
transition  probabilities.  This  problem  has  also  been  investigated  for  the  BPI. 
If  m < 1 and  h^  > 0 It  can  be  inferred  from  results  in  [43]  that 

*^IJ  '^00 

where  is  the  stationary  measure  (which  is  unique  up  to  constant  multi- 

pliers) normalized  so  that  Hq  ~ ^ hence  it  suffices  to  consider  (Pqq^}- 
As  above,  we  are  interested  in  the  case  where  S i>  X^^sitlve.  V/hen  xn  ^ I 
and  < «•  and  some  further  regularity  conditions  hold.  It  has  been  shown 
[29]  that  Pq”^  ^ Kn’*^  (a  = * < m < ••  a full  analysis  was  given  in 

[27].  Some  general  results  and  examples  were  given  for  m < 1 In  [33]  and 
here  we  shall  consider  the  general  results  in  conjunction  with  Conditions  (i)- 


If  m < 1 , Theorem  3 in  [33]  can  be  expressed  as 


pjjg’  H(n) 


where  0 < K < «•  and 


H(x)  s exp  r logh(l -(i)(m^))dx 
*'0 

X 

> exp  r log[l-G(-log  w(m’‘))]dx. 

When  m = 1 , Theorem  4 of  [33]  can  be  written  as  at  (1)  but  with 
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If 

X 

H(x)  * exp  r log[l -G(x)Idx.  (3) 

In  the  following  theorem  we  define  f = 1 when  m = 1 . 

Theorem  15.  Let  m^l  end  (2.4)  hold  if  m®i. 

If: 

(e)  Condition  (1)  holde  then  } iJLiX' 

(b)  Condition  (ii)  holde  then 

Po"o  " 

where  R(  • ) ie  ev  et  infinity; 

^ (c)  Condition  (iii)/(iv)  hold  then 

Pqo^  ~ Kexp(.un*“^X(") 

where  x{*)  le  ev  et  infinity. 

Proof.  For  Conditions  (i)  end  (ii)  it  simply  suffices  to  observe,  using  Lenome 
1 when  m < 1,  thet  multiplying  the  integrends  of  (2)  end  (3)  by  x gives  fane* 
tlons  which  tend  to  -ei^  es  x — we  define  a ® 0 for  Condition  (I).  It 
follows  thet  H(x)  is  'e^-verylng  et  infinity. 

For  Conditions  (lii)/(iv)  it  suffices  to  observe  thet  the  integrend  of 
~logH  (x)  is  -b-verylng  et  Infinity  whence  the  representetion 
H(x)  • exp(-5x^’^XW)  where  x(.)  is  sv  et  infinity. 

Still  essuming  thet  m^  1,  let  A(s)  ® which  converges  for 

j 0 ^ 8 < 1 [43]  end  setisfles  the  equetion 

i A(s)  ■h(s)A(f(s)). 
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Iteration  and  reference  to  (Z.  1)  yields 

For  0 < s < 1 choose  the  Integer  nfs)  for  which  • < ^i+n(s)’ 

n(e)  as  s — ►!-  and 

whence 

A(8)  A- (s-^-l).  (4) 


Thus  once  we  find  the  asymptotic  dependence  of  n(8)  on  »,  we  can  appeal  to 

(l)o(3)  and  /or  Theorem  15  to  obtain  information  on  the  growth  of  A(s)  and 
n 


hence  of  p, 

j*0  J 


Whon  m < 1 we  have  f . 4=1-  <rf(m”^"^)  whence 

nts) 


, n(s)  . , l-fn(s) 

1 - m < Q(8)  < I - m 


(5) 


giving 

n(s)  ~ ^;log(l-8)"^  (e— ►!-). 

Similarly  if  m s 1 we  obtain 

n(s)  < U(s)  < 1 4 n(s)  (6) 

whence 


n(s)  ~ U(8). 

These  results  may  be  applied  to  parts  (a)  and  (b)  of  Theorem  15  to  obtain,  for 
example,  for  part  (b):  If  m < 1 and  Condition  (ll)  holds  then 

A(s)  ~ [ulog(l-8)'^J*^/R(log(l-s)"S 

and  if  m ■ 1 
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A(8)  (U(«))*<^/R(U(8)). 

Since  logH(x)  has  a monotone  derivative  it  follows  that  in  Theorem 
15,  part  (c),  x(*)  integral  representation 

X _ j 

x(x)  = cexp  r e(y)y  <*y 
•'1 

where  0 < c < • and  c{x)“^0  (x— ^•).  When  m = 1 we  see  from  (6)  that 
U(s)  - n(8)  is  bounded  and  some  algebra  shows  that 


(n(s))‘‘^X(U(«))->-0  (s-vl-) 

and  hence  that 

A{s)  ~K”*exp((U(8))‘‘^X(U(«))). 

When  m < 1 similar  manipulations  show  that 

A(8)  ~K"*exp{5[log(l-s)"V‘^X(Clo8<*'»>”^>l 
provided  that  ^l^logCl-s)  - n(s)  remains  bounded  as  s*-^!-.  However,  it 
follows  from  (5)  that  this  is  equivalent  to  the  finitude  of  Q'<1>).  Thus  if 
m<l.  Condition  (iii)  or  (iv)  holds  and  if  E(ZjlogZ|)<»  then  (7)  holds. 


6.  The  reversed  BPI 

Esty  [lO]  derived  some  limit  theorems  for  a time  reversed  simple 
branching  process.  We  now  consider  the  corresponding  reversed  time  BPI. 
For  0 < m < a*  let  h^  > 0 and  R s 1 /h(q).  When  m > 1 , R > 1 and  it  is 
known  that  converges  if 

s < q,  y(q)  * 1 and 

REytPij'yj.  0) 


Thus  if  we  define 


Pj  » yj/yQ  w'i®"  m > I 


then  it  follows  that  for  all  m the 
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strong  ratio  limit  property  holds;  if  i»j>k,  then 

,,  (n-m) , (n)  _m  l-k  . 

n— • 


rjj 

n— • 

1£  0<nj<.  ..<nj^  and  0<lQ,ij Ij^  are  Integers  and  the  l*s  are 


in  S then 


lim  P(X  = I, *n> 

n-n.  1 n-n.  K'  n u 


n-n, 

n 1 


exists  and  is  the  joint  probability 


h \ 

1 k 


of  a Markov  chain  r = TV  1 having  the  n*step  transition  probabilities 

q‘”>  = (2) 

We  shall  call  v the  reversed  BPI. 

VThen  m > 1 the  results  quoted  above  show  that  $ is  poslttve  for  y 
and  the  limiting-stationary  distribution  is  If  m < 1.  R » I and  it 

follows  that  the  classification  of  S ior  y is  exactly  that  for  %•  In  particular 
Theorems  A and  B hold  for  y . We  shall  now  consider  cases  where  j is  not 
positive  for  y . 

It  follows  from  (2.1)  and  (1 ) that  if  0 < t < 1 


SPjt4j(s  “)  » P^"Nt)A(sf^(t)) 

= A(t)A(s£  (t))/A(f  (t)). 
n n 


(3) 


If  1 and  Condition  (i)  holds  or  if  m < 1 and  Conditions  (ii)-(iv)  hold  and 
in  the  cases  of  (iii)  and  (iv)  we  also  assume  that  E(Z|log  Z|)  < « , than  the 
results  of  the  last  section  show  that  a(x)  e A(1  - x S 1)  is  sv  at  Infinity. 
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If  we  choose  a = s In  (3)  so  that  (1 -f  )/(l‘a  ) — > 0,  the  slow  variation 
n n n 

shows  that 

A(s  f (t))/A(f  (t))  ^ A{s^)/A(f  ). 
n n n n n 

Vitali's  theorem  now  yields 

V 

E,(s  ^ A{s^)/A(f  ). 
in  n n 

Now  chooae  s so  that  (1-s  )’ ^ =«’ *(»»(»„))  where  a = (1-f  Clearly 
XI  O 11  n o 

1 < (l-f  )/(l-s  ) and  this  quotient  tends  'o  infinity,  for  if  not 
“ n n 

a((l-«.)'S/aU  )-^l  M n increases  through  some  subsequence  in  K 
n n Y ▼ 

and  this  contradicts  the  definition  of  s . Thus  E, (s  )— ►lA*  applica 

XI  in 

tion  of  Proposition  1 yields 

Theorem  16.  If  m<l,  A(s)  — (a— ►!-)  au|d  A(l-l/x)  is  sv  at  infinity 
then  (lA(l- (1+ V^)'Sl/A(f^)}  has  a limitinR  uniform  distribution  on  (0,1). 

Theorem  17.  If  m x 1 and  Condition  (ii)  holds  then  f(l-f  .)V_1  has  a 

— T 1 . . ■ II 

limiting  gamma  distribution  whose  density  is  (x**^  ^e  *)/r(ai;)  (x  > 0). 

Remark.  If  then  (X^/yn)  and  (V^/ny)  have  the  same  limiting 

distribution. 

Proof.  First  observe  that  A(l-x'S  i«  ay- varying  at  infinity  and,  from 

(3.  3.2)  and  (3.  3.  3)  (1 -f  (t))/(l -f  ) — **l . Now  set  a « exp  (-g  (l-f  ));  it  is 

n n n n 

easily  seen  that  E (exp(-g(l-f^)V^)) —►  (l+g) 

Theorem  It.  Let  m « 1 . Condition  (iti)  or  (iv)  hold  and  a^  » (l-f^)/nG(n). 

Then  a V y. 

n n 
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Proof.  Let  s s exp(-A«  ) and  observe  that  1-s  f (t)  ~ ~ 1-f  and 

n a n n n n 

hence  that  U(8  f (t))  ^ U(f  ).  It  follows  from  (5  1),  (5.3)  and  (5.4),  a mean 
n n n 

value  theorem  and  the  regular  variation  of  G(  • ) that 

log[A(8  f (t))/A(f  )]  ~ G(n)IU(8  f (t))  - U(£  )J 
n n n n n n 

~ G(n)  \J'{t  )(s  f (t)  - f ) 
n n n n 

.w  -eG(n) 

where,  in  the  penultimate  step  we  have  used  the  mean  value  theorem  again. 

The  theorem  now  follows  on  observing  that  (1-s)  ^ U*(8)/U(8)  — >1;  and 

U(f  ) * n. 
n 
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